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Abstract. Wg study some functionals that describe the density of vortex 
lines in superconductors subject to an applied magnetic field, and in Bose- 
Einstein condensates subject to rotational forcing, in quite general domains 
in 3 dimensions. These functionals are derived from more basic models via 
Gamma-convergence, here and in the companion paper [1]. In our main results, 
we use these functionals to obtain descriptions of the critical applied magnetic 
field (for superconductors) and forcing (for Bose-Einstein), above which ground 
states exhibit nontrivial vorticity, as well as a characterization of the vortex 
density in terms of a non local vector-valued generalization of the classical 
obstacle problem. 



1. Introduction 

In this paper we study certain limits of the Ginzburg-Landau model, which 
describes a superconducting object in an external magnetic field, and the Gross- 
Pitaevsky functional, which describes a Bose-Einstein condensate confined in a trap 
and subject to rotational forcing. 

Most prior mathematical work on these sorts of problems has been limited to 
2-dimensional models that are good descriptions, in various regimes, either of 
very flat, thin objects (superconductors or condensates), or of objects that are 
translation-invariant, or very nearly so, in one direction. Important results about 
the 2-dimensional Ginzburg-Landau model, obtained by Sandier and Serfaty in [24] . 
[25] (see also the book [55]) include the characterization of the applied critical mag- 
netic field, below which the ground state of a superconductor expels the magnetic 
field; and above which the superconductor in the ground state is penetrated by 
magnetic vortices; and a description, in terms of an obstacle problem solved by the 
magnetic field, of the limiting density of magnetic vortices above the critical applied 
field. Similar descriptions of 2d Bose-Einstein condensates hold, though they are 
somewhat less well-documented in the literature. 

In this paper we prove analogous results for the full physical problems of 3- 
dimensional superconductors and condensates. In particular, we find a critical 
applied field for superconductors, characterized by the same dichotomy as in the 2d 
case; and we obtain a description, in the supercritical case, of the limiting vortex 
density in terms of a constrained minimization problem solved by the magnetic 
field. This problem is not a classical obstacle problem, but it is a kind of nonlocal, 
vector- valued obstacle problem with an interesting structure. We also establish cor- 
responding results about vortices in Bose-Einstein condensate wave functions, that 
is. ground states of the Gross-Pitacvsky functional with rotational forcing. These 
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results include the determination of a critical rotational velocity, and a characteri- 
zation, in terms of a nonlocal generalization of an obstacle problem, of the limiting 
vortex density for rotations above this critical value. 

We obtain these results from the study of certain functionals whose ground states 
characterize vortex density, and other associated quantities, in limits of sequence of 
minimizers of (suitably scaled) Ginzburg-Landau or Gross-Pitaevsky functionals. 
In the case of Ginzburg-Landau, this limiting functional was derived in a companion 
paper, see U, as a corollary of a general result proved there about the asymptotic 
behavior of a relatively simple model functional. In the case of Gross-Pitaevsky, 
the derivation of the limiting functional, using results of [4], is given in Section [4.1l 

Although we mostly emphasize the analogy between the problems we study here 
and obstacle problems, there are also close connections between our vortex density 
models and total variation models in image processing as introduced by Rudin, 
Oshcr, and Fatemi [23]. (See U for a survey of related mathematical results.) In 
particular, the functional Q derived in Proposition |4l see (|1.22[) . can be viewed as a 
generalization of the Rudin-Osher-Fatcmi model, and in situations with rotational 
symmetry, it reduces to exactly a (weighted) Rudin-Oshcr-Fatemi functional. The 
paper concludes in Section [5] with a discussion of this and some related issues. 

1.1. the Ginzburg-Landau functional. Let be a bounded open subset of M"^. 
A superconducting sample occupying the region Q may be described by a pair (u, A), 
where u is a complex- valued function on fl and ^ is a 1-form on R'^, that encode 
various physical attributes of the superconductor. For example, |up corresponds 
to the density of Cooper pairs of superconducting electrons; dA can be identified 
with the magnetic field; and the superconducting current is given by 

(1.1) Jau := —{udAU ~ udAu), where dAU := du — iAu. 

Stable states of a superconductor in an external magnetic field H^^ex ~ dA^^ex, 
with A^^ex G ^^/ocC^'^' ^^^'^)' correspond to minimizers (or local minimizers) of the 
Ginzburg-Landau functional: 

Jn ^ 4e ^ 

Here the parameter e is related to physical properties of the superconducting sam- 
ple. We will study the limit e ^ 0, with A^.ex scaling so that it will turn out to be 
comparable to the critical value mentioned above. For a discussion of the physical 
relevance of this scaling, see for example (25] ■ 

The model case is a constant external magnetic field, for which we may take 
Ae,ex = ^Ce{xidx^ — X2dx^) for somc real-valued scaling factor c^, corresponding 
to a spatially constant external field H^ ^x = c^dx^ A dx^, which in this example 
points in the 63 direction. 

The functional makes sense for u £ H^{^l;C) and A such that A — A^^^x G 
i/^(R'^; A^M'^). As is well known, the functional is gauge-invariant in the sense 
that for any such (w. A) and for any function (j) such that d(t) S i?f^(R^), the identity 
J"e(u, A) = J"<;(e*'^u, A+d4>) holds. Moreover, (u. A) and {e^'^u, A-\-d(j)) correspond to 
exactly the same physical state, in the sense that all physically observable quantities 
are pointwise equal for the two pairs. 
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Our starting point is the following, which is an immediate consequence of [4], 
Theorem 4. We use the notation 

HliA^W^) := {1-forms A in H\M.^) : d* A = 0} 

which is a Hilbert space with the inner product (A, := /jjg dA ■ dB. We will 

often write H} for short. We also write H^^K^U) to denote the space of p-forms 
on U with coefhcients in the Sobolev space H'^ . 

Proposition 1. Let Vl <ZM.^ be a hounded open set with boundary. Assume that 
Ae^ex G Hl^ci^^R^) "■'^d, that there exists Aex G Hl^^^A'-R^) such that 

(1.3) llJI"^'^^^'^ mHl{A^R^):^{AeH\A^R^):d*A^O}. 

Let (we, Ag) minimize in H^{il; C) x [A^^^x + (A^M.^)] . Then there exists some 
Aq e [Aex + Hi] and vq G L^{Q,; A^M!^) such that dvo is a measure, and such that 

(1.4) 1^^^ -^0^0 weakly in HKA^R^), 

(1.5) 'If' ^1 ^ vo — Aq weakly in L^(A^Q) , 
|loge| 

and 
(1.6) 

tt^Ja^u, —L^d{jA,u + A,) ^ Uvo in W-'^PiA^n) Vp < 3/2. 
|loge| 2|loge| 2 

Moreover, (vo,Ao) minimizes the functional 

(1.7) J^{v,A) = ^\dvm) + ^\\v - Alli.(Aio) + IWdA - H^xWh^A^ 



in L'^{fl; A^M.^) x [A^x + Hi]. Here \dv\ denotes the total variation measure asso- 
ciated with dv. (We understand J^{v, A) to equal +oo if dv is not a measure.) 

Our new results about superconductivity in this paper are derived entirely by 
studying properties of the limiting functional J^; the connection to the more basic 
Ginzburg-Landau model is provided by the above Proposition [T] 

Note that T is gauge-invariant in the sense that if 7 S i/^(A^R^), then 

(1.8) J^{v-dj\n,A + dj)^J^{v,A). 

Our first main result reformulates the problem of minimizing J- through convex 
duality, the relevance of which in these settings was first pointed out in 6J. 

Theorem 2. Let C K."^ be a bounded open set with boundary, and assume that 
Aex e Hl^^{A^R?). A pair {vq,Aq) minimizes T in L'^{A^^) x {A^x + HHA^I^^)] 
if and only if the following two conditions are satisfied: 

1. The 2- form Bq — d(AQ — Aex) belongs to the constraint set 

(1.9) C := e H\A^R^) dH^{AV), supp(d*B) C n, < ^j, 
where 

(1.10) 11^11,:= sup I / B-da : a e H'^ (A^R^) , [ < l|. 
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In addition, 

(1.11) d* + ln{Ao - vq) = m^^, 
and if vq G L'^{K^M?) is any 1-form such that uo|f2 = vo, then 

(1.12) / B^-dvo^-\ f \dvo\. 

2. Bq is the unique minimizer in C of the functional 

(1.13) B ^ £o{B) := I- [ \Bf + ^f\d*B + A,,\^. 

It is clear that if < oo then 

(1.14) / B-dv^O for all u e i7i(Ai]R3) such that in f]. 

Remark also that (|1.12p implies that ||i?o||* = ^ if the vorticity dvo ^ 0. This is 
related to the following necessary and sufficient condition for the vorticity to vanish. 

Theorem 3. Let f2 C ^ he a hounded open set with houndary, and assume 
that € Hl^^{A^R^). Let (uo,^o) minimize T in L'^{K^Vl) x [A^^ + iji(AiR3)], 
and let B^ denote the unique minimizer of £q{-) in the set 

(1.15) C := |b e H^{K^M?) n dH^{k^^^) : B satisfies (1111) 1 

where £q is defined in ()1.13p . 

Then dvo = if and only z/ < i. 

As a result, if (ug, Aj) minimizes T^^ in H^{fl;C) x [A^^^x + -ff* (A"'"R'^)], then 
I \oge\~^JA,u^ as e-^ in W~^'P,p < 3/2, if and only if < i. 

The subtle point in Theorem [3] is the identification of the correct space C in 
which the magnetic field for the vortex-free minimizer is energetically optimal. 

Remark 1. We give a different (but necessarily equivalent) characterization of when 
vq is vortex-free, and a different dual problem, see Theorem [9] and Lemma [10] in 
Section 13.31 

Remark 2. Observe that for B E C we have, by virtue of Hahn-Banach theorem, 

=inf{i|/3iUoo(o), l3eH'iA^R''), d*l3 = d*B}, 

so that can be interpreted as a nonlocal L°° norm oi B G C for O C M.^ with 

iV > 3, while for C we have \\B\\^ = ||B||i=e(o), since in that case in (fll0|) 
one can test with forms aj G i/^(A^R^) with daj = fj{x)dxi A dx2 such that the 
scalar functions fj{x) converge to a Dirac mass 6xo for arbitrarily fixed xq € ^■ 

Hence in the 2-dimensional case C the constrained variational problem 
(|1.9|) , (|1.13p of Theorem [2] corresponds to a classical obstacle problem (as is well- 
known), while in three (or higher) dimensions one may interpret it as a generalized, 
nonlocal, vectorial obstacle problem. Remark also that norms related to || • ||* have 
been studied in the context of critical Sobolev spaces (see [5l [28]). 

Remark 3. Condition p.l4p is easily seen to imply that svipp{d*B) C fJ. The 
converse holds if and only if fl is simply connected. We do not know whether the 
minimizer of £0 in the space C" := {B € H^{A'^R^)ndH\A^R^) : supp{d*B) C n} 
coincides with when fl fails to be simply connected. 



VORTEX DENSITY MODELS 



5 



Remark 4. It may appear unsettling that the functional So contains the non-gauge- 
invariant quantity Aex ■ This can be effectively eliminated, however, by decomposing 
Aex\n = + Al^, where Al^ € (kerd)-^ and Al^ e kerd (see (EIH), (HH)). Then 
Al^ is gauge- invariant, and J^^\d* B + Aex\'^ = J^ld* B + AlJ'^ + \AlJ'^ . 

As mentioned before, a rather complete analysis of the asymptotic behavior of 
the Ginzburg-Landau functional for superconductivity in 2d can be found in [26j . 
In the 3d case, some results in agreement with Proposition [1] have been obtained 
by formal arguments, as in [8], and, rigorously in [3] for the case of the ball, using 
some arguments of [18j . In particular, [3] identifies a candidate for the first critical 
field for the ball in a uniform applied field. Our results show that this candidate 
first critical field is in fact correct. (In fact, it agrees exactly with the alternate 
expression for the critical field alluded to in Remark [TJ see Theorem [S]) Critical 
fields on thin superconducting shells, among other results, have been derived in 
[21 llOj via a reduction to a limiting problem on a 2d manifold. 

1.2. the Gross-Pitaevsky functional. The second main object of study in this 
paper is a variational problem that describes a Bose-Einstein condensate with mass 
TO, confined by a smooth potential a : M'^ — )■ [0, oo) such that 

(1.16) aeC°°(R^), a{x) ^ +00 as\x\ +00, 

and subjected to forcing $e that in general depends on a scaling parameter e. In the 
model case corresponding to rotation about the z-axis, $e := ^Cf{xidx^ — X2dx^), 
and a{x) grows quadratically or faster. In this situation, a stable condensate may 
described by wave function uq that is a minimizer, local or global, in a function 
space to be specified shortly, of a functional that may be written in the form 

(1-17) Geiu) = [ hvu\^-<^,.JU+^{p-\u\Y + ^\u\', 

Jr3 2 4e^ 2e^ 

where j{u) ~ ^{udu — udu). The functions p, w appearing in the functional Q^^ are 

determined by the trapping potential a and the mass m as follows: 

(1.18) 

p{x) := (A — a{x))^ , w{x) := (A — a{x))^ , for A such that / p dx — m. 

The last condition clearly determines A uniquely. 
We study Ge in the function space 

(1.19) Hl{R^;C) := Hi := completion of C;?°(M^; C) with respect to || • |U, 
where the norm || • \\a is defined by := /^a \du\^ + (1 + a)|up. We also define 

Hl„,{W?-C) Hl„, := {u e Hi : j \u\^ = m}. 

We will study the behavior of minimizers of Q^^ in Hi 

Throughout our discussion of the Gross-Pitaevsky functional we will use the 
notation 

(1.20) n = {xem? : p{x) > 0}. 

We will always assume that A is a regular value of a, so that \Da\ > c > on d^l, 
and hence w > in M'^ \ f2, and 

(1.21) |Vp(x)p +p(a:) > c> 0, p(x) > cdist(x,9rj), for all a; € f^. 
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1.2.1. results: Bose-Einstein condensates. Our results for the Gross-Pitaevsky func- 
tional parallel those we obtain for the Ginzburg-Landau functional: we identify a 
limiting variational problem, sec (jl.22p . (|1.23|) below, characterize when minimiz- 
ers of the limiting problem are vortcx-frcc, and obtain a description of minimizers 
of the limiting problem as solutions of a sort of nonlocal vector-valued obstacle 
problem. 

Wc start by proving a theorem that characterizes F-limits of the Gross-Pitaevsky 
functional, see Theorem [T^] in Section 2] . This is parallel to Theorem 4 from [3] 
for the Ginzburg-Landau functional, and the proof relies on results from [?] on the 
reduced GL functional. An immediate consequence of Theorem [T2l is the following. 

Proposition 4. Assume that $c = | loge|$, with $ G Lf„^(A^M^) and that |<f>(a;)P < 
Ca{x) outside some compact set K . 

Assume that minimizes in H\^^. Then 

\u^\ p in L^{R^) 

for p defined in (|1.18l) . and there exists Jq G L^/^(A^fi) such that 

I \og e\~^ju^ jo weakly in L'^^^{M.^) 

Moreover, jo = P^Oj where vq is the unique minimizer of 

(1.22) c;(„):^^p(^M!_„.$ + l|d^;|y 
in the space 

(1.23) LliA'n) |« e Ll,{A'n) : ^ p\vf < ooj . 

(We set G{v) — +00 if dv is not a Radon measure or if p is not \dv\-integrable.) 

Remark 5. We also establish parallel results in the case when $c = ^/gl^ with 
I loge| <C <C e^^. In particular, in this case the limiting energy corresponding to 
([L22| is given by G{v) = p{\v - $p - |$p). See Remark[7]for details. 

It is clear that i^o = $|n is the unique minimizer oi Q. In particular, this implies 
that for a supercritical rotation around the (vertical) X3 axis, corresponding to 

= ^/g^^{xidx2 — X2dxi) with as above, the limiting ground-state vorticity is 
given by dv^ = codxiAdx2, corresponding to an asymptotically uniform distribution 
of vortex lines throughout the condensate, regardless of its geometry or topology. 
This generalizes to 3 dimensions results obtained in |12j in the 2d case. 

We next identify a necessary and sufficient condition on $ and p for minimizers 
of the limiting functional G to be vortex-free, by which we mean that dvo = in 
il. For this result, it is useful to note that the space L^(A"'^ri) defined above is a 
Hilbcrt space with an inner product and norm that wc will write as 

{v,w)p:^ pv-wdx, \\v\\p -.^ {v,vyj'^. 
Jn 

We will sometimes use the same notation to denote the p- weighted inner product 
or norm for /c-forms with values of k other than 1; the meaning should always be 
clear from the context. We let Pp denote the orthogonal projection with respect to 
the Lp inner product, onto (kerd)p, where 

(kerd)p Ll-c\osuie of {(f> € C°°{A^n) ■.d(j) = 0, \\(j)\\p < 00}. 
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We will also write for the eomplementary orthogonal projection. Note that if 
w £ Image(Pp'-) = (kcrrf)^ , then J (pw) ■ cj) ~ for all cj) e (kerii)^ 3 kerd. Thus 
pw £ (ker d)-^ , and so it follows from the standard unweighted Hodge decomposition 
(see section [2?2l and in particular (|2.7p ) that 
(1.24) 



\d*P\ 



2 



Mil- 



V w e (kerd);^, 3/3 e Hl{A'^n) such that w = — and / 

P Jn P 

Thus if $ e L^^ there exists /3$ e iJj^ such that fi*/3$ e ^nd $ = Pp<I> + 
We now state 

Theorem 5. Suppose that Q is a bounded, open subset ofM.^ and that p E C'^{fl) 
and $ e Lfg^(AiM^) n Lp(Aifi) are given, with p satisfying (|1.2ip . 

Lei /3$ G i?^(A^fi) &e sitc/i i/iai Pp'"$ = o,nd let (3q minimize the functional 

1 f 



(1-25) 

m the space 

(1.26) |/? e i7^(A2l7) : ^ELliA'n), - /3$||,, < i 
w/iere 

(1.27) ||/3||p, := sup \ f (3-dw: w£ C°°{A^Cl), f p\dw\ < 1 

Then vq = Pp$ + ^^-^ is the unique minimizer of Q{ ) in Lp(A^fi). 
Moreover, 

(1.28) / _ /3o) . = i /■ p\dv^\. 

Jn ^ Jn 

Finally, dvQ ^ if and only (f < 5. 

Note that (|1.28l) states that the action of the vorticity distribution dvo on the 
potential f3o — is the largest possible given the constraint (jl.26p . Similar consid- 
erations apply to (|1.12p and (|1.9p in the case of superconductivity. 

We finally remark that in general there can exist (3 g Hlf{A^ft) satisfying g 
L'^(A^il) and i|/3j|p* < 00, but such that f3 is not |(ii(;|-integrable for some w £ 
Lp(A^ri) such that / p\dw\ < 1. Hence the restriction to smooth 1-forms w in the 
supremum that appears in the definition of the || • \\p.^, norm. 

1.2.2. related results about 2d BEC. We do not know of any source in the literature 
that establishes 2d results analogous to Proposition |4] and Theorem [S] Such results 
are however in some sense known, at least as folklore, and can be established by 
arguing exactly as in the proofs we supply here in the 3d case, but taking as a 
starting-point results from J9] about F-limits of the reduced Ginzburg-Landau 
functional in 2d, rather than the analogous results about the same problem in 3d 
from [4], which (together with very general convex duality arguments) are the chief 
input in the relevant proofs. 

In particular, limits of sequences of minimizers in 2d are described, in the same 
sense as in Proposition 31 by a functional Q on Lp(r2') of exactly the same form as 
in (|1.22p . for a suitable Q! C M^. More generally, this functional can be obtained as 
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a F-limit of the scaled 2d Gross-Pitaevsky energy, completely parallel to Theorem 
1121 Moreover, this limiting functional admits a dual formulation as functional with 
constraints, parallel to that in Theorem [SJ and from this one can easily determine 
a necessary and sufficient condition for the limiting vorticity to vanish. 

On the other hand, for more extreme rotation regimes in anharmonic trapping 
potentials in 2d, a quite detailed analysis has been carried out recently in [T2 j fTT | [22 ] . 

In a different direction, the critical rotation has been derived in certain higly 
symmetric domains in for example [141 1151 [T] . These references also examine the 
behavior of minimizers for slightly supercritical rotations. 

The main difference between 2 and 3 dimensions is the form of the constraint 
in the limiting variational problem. In particular, in 2d, as in 3d, it is the case 
that if vq minimizes Q, then dvo = d{^^—^), where the potential /3o minimizes the 
functional (|1.25p , subject to the constraint (|1.26p , where the norm in the constraint 
is defined as in (|1.27p . The difference is that in 2d, the potentials /3 are 2-forms 
on M^, and so can be identified with functions. And since it is not hard to check 
that {doj : Jj^ /9|(iw| < 1} is weakly dense in the set of signed measures /i such that 
Jj^ pd\n\ < 1}, the 2d constrained problem reduces to minimizing (|1.25p in the set 

(1.29) |/3 e H'iA^n) : \\1{I3 - /3<,)|U=o < 1| . 

This is a classical (weighted) 2-sided obstacle problem; for many $, using the 
maximum principle it in fact reduces to a one-sided obstacle problem. 

Thus we view the problem in Theorem [5] as a nonlocal, vector- valued analog of 
the classical obstacle problem. 

1.2.3. previous work in 3d. For trapping potentials of the form a{x) — Y^'^i=i '^i-'^i 
with uji > 0, and for — \ loge|$ = A| loge|(a;i(ia;2 — X2dxi), the description of the 
critical rotation given in Theorem [5] was obtained in |17j . building on earlier work 
of [2]. Indeed, these results show that if ||/3<i.|j* < 5, then vorticity vanishes in the 
sense that Jv^ — as e ^ 0. This is stronger than the estimate | loge|~^ J^Je — > 
that follows from Proposition |4] and Theorem [5l The same paper [17] characterized 
the P-limit of the Gross-Pitaevsky functionals for the particular choice of a and <I>e 
described above, in cases where one has an a priori bound on the part of the energy 
associated with the vorticity of the condensate. A similar F-limit result was shown 
by Montero |21| to hold for very general trapping potentials a and forcing terms <I>e. 
Although it was not done in , this result could in principle be used to prove that 
Jve — >■ for subcritical rotations in these situations. The a priori bound on part of 
the energy required for these results means that they cannot give any information 
about minimizers for supercritical rotations. 

Acknovi^ledgments. S.B., R.L.J, and G.O. wish to thank FIM at ETH Ziirich, 
where part of this work was completed, for the warm hospitality. 

2. BACKGROUND AND NOTATION 

2.1. differential forms. If U is an open subset of R", we will use the notation 
W^'P{A''U) to denote the space of maps U — >■ A'^R" (that is, fc-forms on U) that 
belong to the Sobolev space W-^'P. A generic element uj £ W^-'P{A'^U) thus has the 
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form 

(2.1) ujo.dx'^^ A • • • A da;"" 

{a:l<Qi<...<Qfc<n} 

with Wq, S W^'P{^;M.) for every multiindex a. We use the notation Lp{A'^U), 
C°°(A'^C/), and so on in a parahel way. 

For an open set with nonempty boundary and ui € C'^(A'^f2), we define ujt 
and cjTv in C°{A''dn) by 

wt = i*'^, where i : dQ — > f2 is the inclusion; luj^ = uj\gn — i^t- 

One refers to ujt and wat as the tangential and normal parts of w on dQ. We 
will use the same notation wtjWjv to refer to the tangential and normal parts of 
(the trace of) a form w € W^'P{A'^U), which one can define by noting that for 
example the map w i— >■ wt, well-defined on a dense subset of W^'P(A'^ri), extends 
to a bounded linear map W^'^{A''Q) — > LP(A'°i9il), or equivalently by applying the 
pointwise definition of ujt, say, to the trace of w at a.e. point of dQ. 

If are elements of L^(0;A'°R"), written as in (|2.ip . we will write a; • to 
denote the integrable function defined by 

{Q:l<ai<...<Q/^<n} 

This allows us to define an inner product on spaces of differential forms in the 
obvious way. We write d* to denote the formal adjoint of d, so that J duj ■ (j) = 
J Lu ■ d*4> when w is a smooth k — 1-form and 4> a smooth fc-form for some k, and at 
least one of them has compact support. Then 

d*(/) = (— 1)'^ *cJ*(/) if is a fc-form, 

where in M'^, the -k operator, mapping fc-forms to (3 — fc)-forms, is characterized by 

UJ /\ i^(t) ~ i^UJ A (f) ~ UJ ■ 4> dx"^ /\ dx^ A dx^ . 

(In even dimensions one must be more careful about signs.) 
We will use the notation 

W\''P{A^Q) {uj e W^^^{A^n) : ujt = 0}, 

WljP{A''n) := {a; e W^^P{A''n) : ujn = 0}. 

and 

TiTiA'^n) := {uj e W^^P{A''n) -.ujT^O^duj^ d*uj = 0}, 

-HNiA'^n) {uj G W^-P{A''n) ■.ujN^O,duj = d*uj = 0}. 

In fact forms in and "H^ are known to be smooth. 

Gauge-invariance implies that the set of minimizers of J^^ is noncompact in 
-ff^(f2; C) X [Aex + H^{A^M.'^)]. In order to remedy this, we will often restrict to 
a smaller space. Thus we introduce 

(2.2) Hl{W^;A^R^) = Hi := {A e H\W^; A^R^) : d*A = 0} 

with the inner product (A, i3)^i^^i^3j — (A, i?)* (dA, dB) 1^2 (^/^2T^3y This makes 
iJ^(A^]R'^) into a Hilbert space, satisfying in addition the Sobolev inequality 

||^||l6(A1K3) < C'II^IIh1(A1E3) • 
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In view of standard results about the Hodge decomposition, given any any 1-form 
A such that A € A^^ex + H^{M.^; A^M^), we can write B := A - A^^ex e in the 
form 

B^B + d(j), where B e Hi and d0 e ij^(M^; A^R^). 

Thus given any pair (u, A) G H^{il; C)x [A^x+H^ {A^M."^)], there exists an equivalent 
pair (u, A) = (ue^**^, A ~ dcj)) in H^{^; C) x [A^^ex + Hi], so that in restricting J^^ 
to C) X [A^,ex + Hi], we do not sacrifice any generality. 

2.2. Hodge decompositions. We will need several Hodge decompositions. First, 
on a bounded open domain fl with boundary, we have, for every integer k € 
{0, . . . ,n} the decompositions 

(2.3) L^iA'^n) = dH\A''~'^n) ® d* Hl^{A''+'^n) ®HN{A''n) 
and 

(2.4) L'^iA'Tl) = dff|(A'=-if7) ® d*H^{A''+'^n) ® HTiA'^n). 

These are known from work of Morrey (see also [16], Theorem 5.7). The first of 
these, for example, means that every oj £ L^(A'^ri) can be written in the form 
uj = da + d*fi + -f, where a G H^{A^~'^n), (3 e d*7J^(A'=+if]), and 7 G and 
moreover da,d*l3, and 7 are mutually orthogonal in L^. 
We will sometimes use the notation 

(2.5) Vevd = H^{n)®HN{A^^), (kevd)-^ = d*Hl^{Ii^n). 

This is justified by the following considerations. First, we claim that for v G 
L2(Aifi), 

(2.6) dv = as a distribution on fl <=^ v G dH^{n) ® HAr(A^ri). 

Indeed, to prove that v G dH^{U) © HNiA^il), it suffices by (|2.3p to verify that 
w ± d*Hlf{A^n). Fix any G iJ^(A2f]), and let Xe S C;?°(r2) be a sequence of 
functions such that = 1 in {x G : dist(a;,9ri) > e}, HVxelloo < C'e- Then the 
assumption that dv = in f2 implies that 

0= / vd*iXeP)^ f XeV-d*l3+ f V^-idXeA*P) 

Jn Jn Jn 

for every e. Thus 

/ v-d*f3=\im / XeV ■ d* 13 ~ — lim / w • *(dxe A */3) = 
Jn ^^^Jn 

where the last equality follows from the fact that /3n ^ 0. This proves one impli- 
cation in (|2.6p . and the other is obvious. 

Similarly, the Hodge decomposition implies that if w G L'^{A''n), then 

(2.7) UJ e d*Hlf{A''+'^n) /a; •0 = V G i/^(A''r2) with d0 = 0, 

Jn 

We also define 

(2.8) P := L^{A^n) orthogonal projection onto d*W^^{A^fl) = (kcrd)-^. 

Given A G i/^(A^K'^) for example, we will generally abuse notation and write PAi 
instead of P{Ai\n)- We remark that 

(2.9) ll^'S|li2(,,)=inf{||i3 + 7||2 ■.jeH\n),dj^O}. 
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For applications to Bose-Einstein condensates we will need a Hodge decompo- 
sition in the weighted Hilbert space Lp(A''f2). In particular, in the notation from 
the introduction (compare (|1.23|) . (|1.24p ). we may decompos43 </> G L^{A^^) a-s 

(j) = uj + ^, iu^Ppcf>e{keid)p, ^eLliA^n). 

P P 
For refined results assuming p and (j) sufBciently smooth, see |21| . 

2.3. duality. Wc will frequently use the following basic result, whose proof can be 
found for example in }13j . 

Lemma 6. Assume that H is a Hilbert space, and that I : H ^ (—00,00] is a 
convex function and that I{x) < 00 for some x G H. 
Let G{x) ■.= I[x) + \\\x\\']j. 

Let I* denote the Legendre-Fenchel transform of I , so that 

r(e) := sup((e,.T)H-/(x)). 

xeH 

Then if we define (x) I*{—x) + the following hold: 

(1) There exists a unique xq H such that G{xo) = miii/f G'(-). 

(2) The same xq € H is the unique minimizer of G^ in H. 

(3) G{x) + G'^{y) > 0, and G{x) + G{y) = if and only if (x, y) = (.tq, xq). 

3. VORTEX DENSITY IN 3d SUPERCONDUCTORS 

3.1. A dual variational problem. We start with the proof of Theorem [51 in 
which we identify a variational problem dual to that of miminizing T , which de- 
scribes the limiting density of vortex lines in a superconducting material subjected 
to an applied magnetic field. We then use this dual problem to prove Theorem [31 
giving a necessary and sufficient condition for the limiting vorticity to vanish. 

In the next section we present several different and, actually, simpler derivations 
of (an equivalent but different-looking expression for) the critical field. The ap- 
proach presented here, although a little more complicated, has the advantage of 
yielding the dual problem of the statement of Theorcm[2l which clearly generalizes, 
in an interesting way, the obstacle problem identified in the 2d literature, see |26j . 

Proof of Theorem\^ Step 0. Let us write ^ = A\n — v and Q = A — Aex, so that 
in terms of the ^, C variables, 

Hv, A) = \f + \d{C e + + \l MCP F{t C). 

Also, let H :— L^(A^fi) x Note that iJ is a Hilbert space with the norm 

IKCOfff :=lieili2(o) + ll'^CIIi.(E3) 

and the corresponding inner product. We next define 

■.= lJjd{c-^ + A,,)\ 



Notice that our notation is inconsistent, with Pp = projection onto (kerd)p for Bose-Einstein, 
and P := projection onto (kerd)-'- for superconductivity. These conventions are convenient how- 
ever, and we do not think they can lead to any confusion. 
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SO that 



As usual, / is understood to equal +00 if d{C, — C + A^x) is not a Radon measure. 
Let /* denote the Legendre-Fenchel transform of /, so that 

r(c,c)= sup {((c,c),(r,c))H-/r,c)} 

Let us further write 

step 1. As remarked in Lemma [S] above, (CoiCo) minimizes F if and only if 
(— ^0, ^Co) minimizes . To compute /*, note that for (^, C,) G H, 

r(^,C) = sup |((^,C),((r-^e.)+Ae„C*))H-^ / M(C-(r-Ae.))| 

(3.1) 

= (e,Ae.)L2(n)+ sup \{{^X),iC,C))H-l [ \diC ~e 

It is clear the supremum on the right-hand side equals zero if (^, C) satisfies 

(3.2) / l^^-C+dC-dC <l [ \diC-e)\ forall(r,C)ei?. 

and if this condition fails to hold, then (by homogeneity) the sup in p.ip is infinite. 
Thus 

r*(/^ /-) = f (C,^ea;)L2(a) if holds 

\ +00 if not. 

It follows that 

ptr^n = / ^ll(^ + ^-'Olllf-^Pex|li.(o) if (E3 holds 

^^'^^ \ +00 if not. 

Step 2. We want to rewrite F^ in a more useful form. To this end, we first 
claim that (^, C) € H satisfies p.2p if and only if 

(3.3) / dC-dC < If \dC\ for a\\ C & H\AV) 
and 

(3.4) C e H'L n Hi and d*dC + loC = 0. 

Step 2a. First assume that (|3.2p holds. Note that since (^, Q ^ H 

{i^X),{CXl)H<l f \diC~C)\ foraU(^,C)ei'(Air!)xi7i(AiR3). 



This follows from dXH, since we can write (r,C*) e i^(f2) x H'^{R^) as (r,C') + 
(0, C") with (r, C) G and C" -L Hi, so that dC = 0. 

Now we immediately obtain p.3p by taking (^*,C*) of the form (0,C*) in the 
above inequality. Similarly, by choosing (^*,C*) of the form ±((^*|a,^*) we find 
that 

(3.5) / {dC-dC + ln^-C) = for all C* e fl"^(AiM3). 
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Since d*C = for all C e Hi, we see from p.5|) that — A^ + lsiC = as distributions, 
and hence from elliptic regularity that C € i/;Q^(R'^) and that d*d(^ + ln£, ~ a.e. 
in R^, so that holds. 

Step 2b. Conversely, suppose that ((331) . ([Q]) hold. Clearly implies ((33|) . 
so for (X*,C*) G iji(AiM3) X iji(AiR3), 



U-x* + [ dC-dC^ f dc- d{C - X*) 

t33t 1 /■ 

< - jjd{C~x*)\. 



Thus p.2p follows whenever ^* is the restriction to D. of some X* £ i7^(A4 

We next deduce from this that p.2p holds whenever ^* g L^(A^r2). We may 
assume that d^* is a measure, as otherwise the right-hand side of p.2p is infinite 
and there is nothing to prove. Then, given (^*,C*): it suffices to find {X*,Q) S 
H\A^R^) X iji(Ai(R3) such that 

X*\n C weakly in L'^{A'^n), 

^2 dC dC weakly in L^(A^R^), and 



id(c-x:)i- / M(c-r)i. 

Jo 

To do this, we start by fixing, for all e sufficiently small, a diffeomorphism 
: M''' ^ R3 such that 

(3.7) *,({.T e R3 : dist(x, fl) < e}) C fl, = a; if dist(x, 90) > 
and 

(3.8) |li5«'.-/||oo < CV^, \\D^7' -I\\oo<CVe. 

For example we may take in {x € R'^ : dist(a;,90) < y/e} to have the form 
^e{s) = X — fe{d{x))D{x), where P{x) is the outer unit normal to 9fi at the point of 
dn closest to X, and d{x) is the signed distance (positive outside fi, negative in n) 
to dn, and /e is a nonnegative function with compact support in (— -v/e, -v/e) such 
that |/,'| < Cy~e and /,(e) > e. 

Next, let ^* denote some extension of ^ to an element of L^(A^R'^), and let ip^ be 
a smooth nonnegative radially symmetric mollifier with support in i3(0,e/2) and 
such that / (j>e — I- 

Then we define 

X: :=i:,*{^tC), and C := * (*f C)- 
Then the verification of (|3.6p follows by a reasonably straightforward, classical 
argument. (See for example the proof of Lemma [T3l at the end of Section l4Tl where 
similar computations are carried out in detail in a somewhat more complicated 
setting.) 

Step 3. In view of p.4p . we can eliminate ^ from the expression for to find 
that 

FUf^n-i lllM*^^C + Ae.,C)lll^-iPe.|li2(^) if (EH) hold 

\ +00 if not. 

We now rewrite everything in terms of ^ = (^+Aex, v = ^|o— ^ = (C+^ea;)|n — 
and B = d{A - Aex) ^ dQ. 
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First, the constraints p.3p . p.4p are equivalent to the conditions appearing under 
part 1 of the statement of the theorem, that is, 

(3.9) BeC, and d* B + ln{A ~ v) ^ 0, 

where the constraint set C is defined in ()1.9|) . 
Second, it fohows from Lemma [5] that 

{vo,Ao) minimizes J" in L'^{Q) x [Aex + V-o] 

{Aq — vq, Aq — Aex) minimizes F in H 

{vq — Aq, Aex ^ Aq) minimizes in H , 

■4=^ (wo, Aq) minimizes in L^(17) x [A^x + "Ho], 

where F\v,A) := F^{v - A,Aex - A) + ^AexWl-n^^y 
Thus 

(310) j-t(„ ^) ^ / + if dHI) holds, and 

Rewriting 7-"^^ in terms of -B := d{A — A^x), it follows that {vq,Ao) minimizes F^ 
if and only if conditions 1 and 2 from the statement of the theorem hold. Remark 
finally that (|1.12|) follows from ()l.lip and the relation 

(3.11) i / \dvo\ + [ {vo-Aq)-vo^O 

which in turn follows by stationarity of F{vo,Ao) with respect to variations Vt = 
e* • Wo around t = 0. □ 

3.2. first characterization of the critical applied magnetic field. We next 
want to prove Theorem [3l which gives a necessary and sufhcicnt condition for the 
vorticity of a minimizer of F to be nonzero. Recall that this criterion involves the 
minimizer of an energy So in a space C, defined in (jl.lSp . We first establish some 
facts about C. Given a function v defined on Q, we use the notation Iqv to denote 
its extension to the function, defined on M'^, that vanishes away from fl. 

Lemma 7. Assume that B e H\A'^E.^) D dH^{A^M.^). 

IfBe C , then supp(rf*B) C and {d* B)\n e d*Hl;{A^n) = (kerd)-^. 
Conversely, given any (j) G d* Hlf{A^fl), there exists B^ G C such that d*B^ — 

Finally, 

(3.12) ^' ^ f] ^^'^ ^ n 

1<P<2 |<9<6 



Remark 6. The proof will show that B^ = d{—A) ^(1^0), where (—A) ^ de- 
notes convolution with the fundamental solution for the Laplacian on M^, with 

(-A)-i(in0) e ni<p<2W^2'P c nr>3Lr 

Proof. Step 1. We first claim that if S e C , then 

(3.13) / d*B-v^O for aU e L2(-yYiK3) g^^j^ ^jj^^ ^ in O. 

(Recall that by definition of C, this identity holds for v E H^{A^M.^) such that 
dv ^0 inn.) To see this, define a diffeomorphism : as in (P?7)l . (j^ . 
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i)= f v,-d*cj)= [ 


hi 





and let denote a symmetric approximate identity supported in -0(0, e/2). Given 
V e L'^{AV) such that dv ^ in n, define Ve := ipe * (^'f w) e H^{A^R^). Clearly 
Ve V in L^(M'^), as e — > 0, and we also claim that = in f2. To see this, note 
that for any (j) € C^iA^n), 

'i>fv{ipe*d*(p)^ / *f w A (-0e * 0). 
Since ^-f w A (^/j^ * 0) = A (*7^)#(d* (V'e * 0))], it follows that 
(3.14) [ dve-cp^ [ V A d{'^-^)*{k{iP^ *(!))).= I vd*-k(j), 

for := (^'7^)^(*(?/'c*0)). The definitions of ^'e and i/'e imply that cpe has compact 
support in fi. Thus /^^ • = J^dv ■ (pe = for every -0 S C^(A^r2), and it follows 
that dve = in 17. Then if B G C, 

d*B • i; = lini / B- dv^ = 0. 
•^^0 Jr3 

Step 2. Now for B G C, if x S C^°°(R3 \ O), then xrf*^ e L^{A'^R^) and 
d{xd*B) in r2, so /jgaxM*^^ = by (I^TO)) . Hence supp(d*B) C Q. Then 
p.l3p implies that for any v G L^(A^f7) such that = 0, 



= 



/ d*B-{lnv)= [ d*B-v 
JR3 Jn 



Thus {d*B)\ii G (kerd)-L = d*Hlj{A^n). 

Step 3. Now, given G d*H]^{A^Lo), let i/- = (-A)-i(lf20), and let := dip. 

Then the fact that cp € d*H'^ implies that d*(ln0) = on R^. It follows that 
d*tp = 0, and hence that d* Bcf, ~ d*dtp = {d* d+dd*)i/j = —Ai/j = ln4'- In particular 
supp(d*V') C ^. 

Finally, to see that B^ G C, observe that Incp G ni<p<2i-''(A^R'^), so elliptic 
regularity and embedding theorems imply B,p G ni<p<2VF^'^ C n3/2<q<6-^', -S^ S 
i/i(A2M3)ndiji(AiM3)^ and in addition ([332]) holds. It is clear that supp(d*B0) = 
supp(ln0) C f2, so B^ G C 

□ 

We need one more easy fact about C . 
Lemma 8. If Bi,B2 G C , then there exists tpi G ni<p<2 

M^2,P(AiR3) such that 

dtjji = Bi and 



[ B,B2 = f A-d*B2. 
Jr3 Jn 



Proof. Let ipi = (—A) ^d*B, so that in view of Remark[6l 

Bi-B2 = / dV^i-S2 = / yJi-d*B2 - / yJi-d*B2 



n 



where the integration by parts is easily justified in view of the decay properties 
recorded in p.l2p and Remark [S] □ 

Now we give the 
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proof of Theorem[^ Step 1. We first assume that < ^- Then, recalhng 

(jl.lip . and recalhng that dBo = 0, we must show that 

= dvo = dd*Bo + dAo = dd* Bq + Bq + Hex in ^■ 

Since Bq and B^ minimize So in C and C respectively, and since C C C, it is clear 
that B^, = Bq if and only if e C, which holds if and only if ||-B*||* < 5. So it 
sufhces to check that 

(3.15) dd*B., +B.,+Hex=0 in n. 

To do this, we take first variations of £0 in C' to find that B^ satisfies 

(3.16) / B^ ■ B + I {d*B^ + Aex) ■ d*B ^0 for all B e C 

By Lemma |8l we may rewrite this as 

/ (Vj + d*B^ + Aex) ■ d*B = for all B e C. 

Jn 

where ip = {—A)^^d*B^,, so that dip ^ B^,. Then we conclude from Lemma [7] that 
{t/j + d*B* + Aex)\n e ((kerd)-L)-L = kerd, and hence that d(V' + d*B^ + Aex) = 
in n, which is ([XTS]) . 

Step 2. Now we assume that dvo = in fi. We will show that in this case, 
£o{Bq) = fo(-B*)- Since B^ is the unique minimizer of Sq in C and Bq £ C C C, 
this implies that Bq = B^, and hence that B^ £ C. 

First note that 

(3.17) / vo-d*B^O for any B e C, 

since d*B e d* H'j^{A'^fl) = (kerd)-'- by Lemma [T] Applying this to B ^ Bq and 
recalling that vq = d* Bq + Aq ~ d* Bq + (Aq — A^x) + Aex in 51, we obtain 

= / {d*Bo + [Ao - Aex) + Aex) ' d* Bo 
Jn 

= / {d* Bo + {Ao ~ Aex) + Aex) ■ d* Bo 

= / \d*Bo\^ + \Bo\^ + Aex-d*Bo. 

(The integration by parts is easily justified using p. 121) .) Using this to rewrite the 
definition of £0 yields 

(3.18) £o[Bo) f Aex- d*Bo + |Ae,p 

^ Jn 

Step 3. Next, taking B* as a test function in p.l6p . we obtain 

/ |B,|2 + lj,(|d*B,|2+Ae, •d*S) =0. 

It follows that 

(3.19) £o{B.) - \ [ {d*B,-Aex + \Aex\^) 

^ Jn 
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3. From p.l6p we also have 

(3.20) 0=/ B, -30 + ln{d*B,+A,^) -(1*80^0 
On the other hand, again using p.l7p . we compute 

= [ vo- d*B, = [ {d*Bo + (Ao - A,,) + A,,) ■ d* B, 

= / (d* Bo + [Ao - A^^) + InA^^) ■ d* B, 

= / d* Bo ■ d* B, + Bo ■ B, + InAe:, ■ d* B^ 

recaUing that d{Ao — Aex) — Bo- And by comparing this and p.20p . we find that 

/ A,,. ■ d*B, = / Ae, • d*Bo. 

This, together with p.lSp and p.l9p . shows that 8o{B^,) = £o{Bo), completing the 
proof. 

□ 

3.3. an alternate characterization of the critical applied field. Our next 
result gives a different characterization of the critical field. 

Theorem 9. Let {vo,Aa) minimize T in L'^{^\h^M?') x {A^.^ 
Further, let 

(3.21) £^IA) = 1 f ln\PA\^ + \dA-H,x\^dx, 

where P is defined in (j2.8p . and let Ai minimize £i in A^x + i?"'^(IR^; A-'^R'^). Let 
ai e H'^{n;A'^R^) be such that 

(3.22) d*ai = PAi, dai =0 in fl, cti,N = on dil. 

(Such an ai exists by definition of P.) Note that Ai and hence a\ depend on A^x- 
Then dvQ ~ if and only if 

(3.23) ||ai||„ := sup dvai<l/2. 

\dv\{n)<i Jn 

Moreover, if dvo ~ then Ao = Ai. 

3.3.1. Theorem\^ via a splitting of T . We will give three proofs of this theorem. 
We first present the most direct proof, which does not use convex duality at all. 

First proof of Theorem\^ RecaU from ((2^ that for v G L^(A^ri), = in if 
and only if w £ dH^{n) ® HAr(Aif7) = kerd, see Define 

f{v, A) := inf { ^"(1; + 7, A) : 7 e ker d} 

/ \dv\ + \P{v- A)\^dx + / \dA-Hex\^dx . 
Jn Jr3 

(Note that the definition p.2ip of £1 can be rewritten £i{A) = 7^(0, A).) It is clear 
that 

(vq, Ao) minimizes J-' <==^ {vq + 7, Ao) minimizes J- for some 7 G kerd. 
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Since we are interested here in dvo, we may consider J-' instead of J^. We rewrite 

J'{v,A)^£i{A) + l I \dv\ + \Pv\^ ~2Pv ■ PAdx. 
2 Jn 

Since Ai minimizes £i, 

(3.24) / InPAi ■ PB + {dAi - He^) ■ dB dx = 
for all B e H\R^;W), so that 

£i{Ai +B)^ £,{Ai) + \ j ln\PB\^ + \dB\'' dx 
for B as above. Given any A, let us write A = Ai + B. Then 
:F{v,Ai+B)=£i{Ai) + l- I \PB\^dx + [ \dB\^ dx 

+ - I \dv\ + \Pv\^ - 2Pv ■ {PAi + PB) dx 
2 Jo 

= 1 / \P{B - + 11 \dB\' + I \\dv\ - Pv ■ PA,. 
^ Jn ^ Jk3 Jn ^ 

For a\ as in the statement of the theorem, 

(3.25) / Pv ■ PAi = Pv d*ai dx = dv ■ a, dx 
Jn Jn Jn 

where the boundary terms arising from integration by parts have vanished due to 

the fact that ai^w = 0. Thus 

(3.26) f{v,A)^£,{A,) + l f \dB\^ + ln\P{v-B)\^+ [ {l-\dv\-dva,). 

^ JR3 Jn ^ 

If condition (|3.23p holds, then J^^{^\dv\ ~ dv ■ ai) > for aU v e L'^{^), and thus 
■F{v,A) > £i{Ai) for all {v,A). Moreover, if {vo,Ao) = {vo,Ai + Bq) attains this 
minimum, then i |(ii?oP + ln|-P(wo — ^o)P = 0, and this implies that dvo = 0. 

And if p.23p fails, then there exists some vi such that —dvi ■ a\ + < 0, 

and then it is clear that T{\v\t Ai) < £i{Ai) = T{0,Ai) for all sufficiently small 
A > 0. Thus J'{vq,Ao) < £i{Ai) for any minimizing {vo,Ao), and then ()3.26p 
implies that dvQ ^ 0. 

Finally, it dvQ = then it is clear from p.26p that J^{Q, Ai) ~ £i{Ai) = min J^, 
and and hence that Aq = Ai. □ 

3.3.2. Theorem\^ via partial convex duality. We next prove Theorem[9]by a duality 
computation that differs slightly from the one used in the proof of Theorem [2l The 
result of this computation is summarized in the following 

Lemma 10. Let 

(3.27) N:=I^CeL\n):iC,OLHn)<\Jm for all ^ e L'' (n}Y 
and define 

(3.28) THA) i Jj\d{A ~ A,,)\^ + ln\A\^) dx ~ ^ disti.(,,)(A, N) 
Then (vq, Aq) minimizes T in i^(fi; A-'^M'^) x \Ai^,j, + HW if and only if 
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1. Aq minimizes in [Aex + Hi], and 

2. ^o|f2 -v^e N, and \\Aq - i'o||L2(n) = distL2(n)(Ao, A^). 

It is clear from the definition that N C (kerd)-'- = Image(P), and it follows that 
^\A) = \ [ i\d{A~A,,)\' + ln\PA\^) dx-idisti. (PA,7V) 



(3.29) 



1 - .2 



= £iiA)--dmti,^a)iPA,N). 



Proof. Wc will compute the convex dual of J- with respect to the "v" variable only, 
treating A as a parameter. Thus, let ^ — Ajo — v, and write F{^; A) — J-{v, A), so 
that 

^ Jn ^ JR3 

Let 

iitA):^l f \diA-0\+CA. Ca:=\I \d{A-A,x)\\ 

ThenF(e;A)-/(e;A) + i||?||i. 
Next, let 

~r{C;A) := sup f(r,OL^(o) -/(C;^)}• 
A short computation like that in the proof of Theorem [2] shows that 



(3.30) i*{C;A) 



+00 if not 



for TV as defined in ([HTTf)) . Now let 

F\C,A) :^r(-r;^) + ^||eir 



(-r,A) + i|ieip-c^ ifreiv 

+00 if not 

hW^* - A\? - h kMA~ A,x)\^ ^ InW) dx -liC eN 
+00 if not. 



Then it is clear that 



-ipFHC:A)=THA) 



as defined above, and that the infimum is attained by a unique ^, the closest point 
to A in the (closed convex) set N . Recall from Lemma[6]that ^ minimizes F{-\ A) if 
and only if it minimizes F^{-', A), and moreover that mine F{-] A) = — min^ F^{'] A). 
Thus 

minJ^(u;y4.) = minmin J^(t;; A) = min(— min/'(^, A)) ~ minJ^"''(A), 

A.v A V A ^ A 

and (uo, Aq) minimizes T if and only if Aq minimizes J-^ and = ^o|o ~ ^Oi where 
^0 is the closest point in N to A. 

□ 

Now we use Lemma [TU] to give a 
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second proof of Theorem\^ Fix (wo,^o) minimizing and Ai minimizing £q. By 
Lemma fTOl Aq minimizes . 

We write — Aq\q_ — vq as above, so that is the closest point to Aq in N. 
From the definition (|2.8p of P we know that cIPAq = dAo, so that 

dvo=0 ^ d{PAo - Co) = in V'{n). 

Also, PAa - Co e (kcrfi)-L = d*M^^'^(A2f7), since the definitions of N and P imply 
that TV C (kerd)-^ and Image(P) = (kerd)-^. Then ([2^]) implies that d(PAo-Co) = 
in V if and only if PAq - Co e d* VF^'^(A2l7) n {dH\n) © ?^(Air2)) = {0}. In 
other words, dvo = if and only if PAq = Co- But since Co is the closest point in 
N to Aq, and hence to PAq, we conclude that 

(3.31) dvQ = PAa e N. 

Next, note that ^ i-> ^^■'■(Aex+A) is strictly convex in Hi, so that the minimizers 
Aq of J"* and Ai of £i are unique. Also, p.29p implies that if PA e N and A' e 
then 

+ A') <£^iA + A') < fHA) + C\\A'\\l.^,,^ < ^HA) + C\\A'\\%^. 

Thus any critical point A of £i such that PA e N must also be a critical point of 
J^^, and conversely. 

(3.32) PAo e if and only if PA^ G A^. 

It follows along the same lines that if dvQ — then A^ ~ Ai. Finally, recalling the 
definitions p.22p of ai and (|3.27p of A^, and integrating by parts as in (|3.25p . we 
conclude that 

PAieN ^ ||ailU<^. 
By combining this with p.31|) and p.32p . we conclude this proof of Theorem |9l □ 

3.3.3. Equivalence of Theorem\^ and Theorem\^ In Theorcm[3]and Thcorcm[9l we 
have derived two necessarily equivalent but rather different-looking necessary and 
sufficient conditions for the vorticity dv^ of a minimizing pair {vq, Aq) to vanish. In 
this section we elucidate the connection between the auxiliary functions P* , defined 
in Theorem [3l and ai, defined in (|3.23p . 

Proposition 11. If dvo = in then PAi — d*ai — ((i*P*)|n, and ||q;i||** = 
||P*||*. 

This can be seen as a third proof of Theorem [HI 

Proof. If dvQ = 0, then Pvq ~ 0, and we have seen that = and Bq = B^,. As 
a result, 

(d*P.)|n = {d*Bo)\n = P{d*Bo) = P{Ao - vo) = PAq = PA, - d*a, 
by Theoremll LemmaEl and pTTT]) . Thus for every v G ^^(A^M^)^ 

/ dv ■ B^, = / V ■ d* B ^ / V ■ d* B = / v ■ d*ai ~ / dv ■ ai. 
Jr3 Jr3 Jn Ju Jn 

Now the conclusion follows from the definitions of the norms || • j|* and || • ||**, see 
pTTOl and ({3:231) . □ 
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j = pv weakly in i'*/^(R^), and v € ip(A^17) with / p\dv\ < oo 



4. VORTEX DENSITY IN 3d BOSE-EINSTEIN CONDENSATES 

In this section wc use the results of [4] to prove convergence as e ^ of Gross- 
Pitaevsky functional Ge, defined in (|1.17p to the limiting Q, defined in (|1.22p . We 
also establish some results describing minimizers of Q. 

4.1. F- convergence. Our first theorem makes precise the sense in which C/ is a lim- 
iting functional associated to the sequence of functionals {Qe)e£{o,i]- The statement 
of the result uses some notation that is introduced in Section [1.2l 

Theorem 12. Assume that = |loge|<I> for $ e Lf^^.{^'^R-^) and that < 
C{a{x) + 1) for all X eR^. 

(i) compactness; Assume that (Me)eg(o.i] C Ham '^^'^ ^^"^^ there exists some 
C > such that 

(4.1) GeM <C\loge\'^ for all ee {0,1/2] 

Then there exists j £ L'*/'^(A"'^R'^), supported in Q, such that if we define v — j;j\n 

and pass to a subsequence if necessary, we have 

(4.2) 

|loge| 

(ii) : lower bound inequality; There exists a sequence of numbers (k^) such that 
if we assume above hypotheses and (|4.2[) . then 

(4.3) liminf | logej-^ {g,{ue) - n,) > g{v). 

(iii) upper bound inequality Given any v £ Lp(A^fi) such that dv is a measure 
on n with p\dv\ < oo, there exists a sequence £ ^ such that (j4.2p holds 
and lime_>o I loge|"^ (5e(ue) - ^e) = Giv). 

The theorem states that the functionals \ \oge\~'^{Ge{-) — Ke) converge to G in 
the sense of F-convergence, with respect to the convergence (|4.2p . As remarked 
in the introduction, Proposition |4] is a direct corollary of Theorem [12] and basic 
properties of F-convergence. 

Remark 7. In fact we prove a more general result than Theorem [121 since we also 
allow higher rotations = VffT^' '^ith I logep ^ ffe ^ e~^. In fact we show that 
for such if Ge{Uf) < Cg^, then after passing to a subsequence, j weakly 

in L'*/3(AiR^), with j = Inpv for some v £ L'^p{A'^Q), and 

(4.4) g-' iG,{-) - «,) A Gi-) , where G{v) = p (^y - $ • z;) . 

As noted in Remark[S]in the Introduction, it easily follows that for rotations around 
the z axis of order | loge| <^ ^Jgl <C e~-^, ground states exhibit an asymptotically 
uniform distribution of vertical vortex lines, generalizing 2d results of |12) . 

The proofs rely at certain points on Theorem 2 in [4]. 

Proof of Theorem[TE and Remark\2^ Let Ge{ue) < Cg^, for |logep < 5e <C e^^, 
and let i>£ = Jgl^- 
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Step 1. First we control the potentiahy negative term in Geiue). To do this, 
recall our assumption that |$P<C(a + l). Since \ju\ < \u\ \du\, it follows that 

\^.-ju,\ < lgemM' + l\du\^ < Cg4a + l)\u\^ + ^\duf. 
But p.lSp implies that a{x) = w{x) — p{x) + A < w{x) + A, so it follows that 

• ju,\ < Cg,{w + A + + ^\du,\\ 

Integrating this over K'^ and recalling that ||it£||2 = rn, we obtain 

/ |$,-ju,| < Cg,K) + Cg,(A + l) / \u,\^<Cg,. 
It follows from this that 

(4.5) / l\du,\^ + l^(p-\u,\^f + ^\u,\^=g,iu,)+[ $,.j«,<C5.. 
Jr^ 2 4e^ 2e^ J^3 

In particular p in L'^(R'^). 

Step 2. Next, 

ll«e||i4(R3) = II Iu^Yl^ < C{\\ \u4' - p\\l. + llplli.) <C + e^GM < C, 

and Step 1 implies that ||(iMe||i2 < C^fgl. Since ||jitc||^4/3 < ||Me||i4||(iU(;||£,2, we 
conclude that {-^juj} is uniformly bounded in _L^/'^(R^), and it follows that there 

exists some j G L'^l'^{^) such that 

(4.6) j weakly in L'^^^ along some subsequence. 

Step 3. Now let denote the minimizer in Hl{R^; M) of ^e(-), where -ffi(M^; M) 
is defined by analogy with H^{M.^;C), see (fTTOl) . 

Note that when / is real-valued, jf = 0, so the forcing term $ • j f vanishes on 
_ff^(M^;M). It is standard that /e does not vanish, and we will assume that /e > 0. 
Then for any u G we may define U := and it is known that 

(4.7) g,{u) = GM.U) = GMe) + H,{U) + / . j[/ , 
where H^{U) = Hg{U;M.^), and for a measurable subset A C we write 

(4.8) H,iU; A) := ^ ^\dU\' + ^{\U\' - if dx. 

See [H] for a proof in exactljQ the situation we consider here, following ideas that 
originated in [50] and have been used extensively in the literature on Bose-Einstein 
condensates. It is also known that — )- p uniformly in R'^, see again |21j . 

Step 4. Now let il' denote an subset of f2 such that il' CC fl, so that p > 2c' 
in fl' for some c', and hence > c' for all sufficiently small e. Let Ue = u^/ fe, and 
note that Step 1 implies that He{Ue) < Cge- Thus the functional 

He{U,)^ f ]-\dU,\^ + \w{u) < C'g,, where I^(C/) = ^(|C/|2 - 1)2, 

verifies in f2' hypothesis {Hq) of Theorem 2 in [4], for q = A. Hence Theorem 
2 of [4] (or arguments such as those in Steps 1 and 2 above) imply that there 

^The Gross-Pitacvsky integrand is written in a slightly different way in in 1211 . but this is 
purely a cosmetic difference. 
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exists some v' in L^(A^ri') such that, after passing to a further subsequence if 
necessary, v' weakly in L'^/^{K^Vl'). Since = fJJ^, one easily checks 

that jUf^ = jUe, and then it follows from (|4.6p and the uniform convergence 
— > p that j = pv' in V,' , and hence that v' = j/p =: w in ^2', independent of V,' . 
It also follows that the chosen subsequence is independent of Jl'. Let moreover 

Me J QMC/ep + liy(u)^ dx, 

be the energy density of iff, and notice also that J^, p^ is uniformly bounded, so 
that after passing to a subsequence, we may assume that there exist a measure 
fi' such that pe po weakly as measures in ft' . It then follows from Theorem 2 
and Remark 4 in [i], that, in the case < C\ logep, po > dx + \dv\), in the 

sense that \dv\ is a Radon measure, and 

Po{U)>lJ^{\v\^dx+\dv\) 

for every open U C ft', while for |logep ^ (/e ^ we have po > ^\v\'^ dx. In 
either case we deduce (using basic facts about weak convergence of measures) that 

(4.9) liminf —i7,(C/e; 17') = liminf / /^Ve > / PMo , 
which yields 

(4.10) liminf —iJ,(C/,; 17') > i / p(\v\'^ dx + \dv\) if .g, < CI log eP, 

9e 2 Jjj, 

and 

(4.11) liminf -i?,(C/,; 17') > i / pli^^dx if | log < < e-2_ 

^-^0 9e 2 

Since this holds for all 17' C 17, it follows in particular that v € i^(A^17) and (in 
case ge < Cjlogep) that dv is a measure on all of 17 with J^^p\dv\ < oo, nearly 
completing the proof of (|4.2p . (We still need however to prove that j is supported 
in 17.) 

Step 5. We next claim that 

(4.12) - I .jT/, -> / p^-v. 

To prove (|4.12p . since $e = ^/g^^ and -^ju^ = -^fe jU^ ^ j = pv weakly in 
L'*/3(17), it is clear that 

- / f^<^>,.jU,^ f p<i>-v, 
and we only need to show that 

Since p = outside 17, in this set we have 

< K\ \du,\ < + ^\du,\^/' = - pf + 
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whence, for any compact X C M'^, we see from (|4.5p that 



f \ju,\<c(e^_ + {e^,y 
V9e JK\n ^ 



^1/3 

lK\ 

Thus -> in L^{K^{K\Vl)) for any compact K. This imphes that j = outside 

il, so that the identity j = pv holds in all of R'^, finally completing the proof of 
()4.2p . and it also implies that 

1 /" r iUf 

' <i>, ju^^ $ • — ^ 



9e JK\n J K\n 

for A' compact. Next, due to (|1.16p . (|1.18p and the assumption that |$p < C(a+1), 
we can find a compact K such that |$P < Cw outside of K , so that (arguing as in 
Step 1) 

1$ • jue\ < + — |up outside of K. 

It follows from this and Step 1 that 

- / \'^e-jUe\ < Ce. 

By combining these inequalities, we obtain the claim ()4.12p . 

Step 6. We now complete the proof of the lower bound inequality. Note that, by 
combining (|4.12p with (|4.10p and recalling (|4.7p . we find that, in case 5e < C| log ep, 

for any open fl' compactly contained in fl. Taking the supremum over all such 51', 
we obtain ()4.3p with = Ge{fe)- Analogously, in case |logep <C <C e^^, using 
(|4.1ip in place of ()4.10p we obtain the lower bound part in (|4.4p 



hminf 1(^,(7.,) -^,(/,)) > / p(\^-<i,.v 

Step 7. Let us prove the upper bound inequality in case < C\ logep. The 
proof in the case | logep <C ^ e^^ follows the same lines and hence is omitted. 
We will use the following lemma. 

Lemma 13. Suppose that v G L^{A^Q) and that dv is a locally finite measure with 
p\dv\ < oo. Then for every S > 0, there exists vs S {A^M.^) such that 

(4.13) [ p\y^^y\^ < J, I p\dvs\ < f p\dv\ + 6 



The proof is given at the end of this section. Now we use the lemma to complete 
the proof of the theorem. 

Fix V and vs as in the statement of the Lemma. 

It is proved in [21j . Lemma B.l, that the (positive) function appearing in 
the decomposition (|4.7p satisfies ||/e||Loc(]j3) < cq, for cq independent of e € (0, 1], 
and moreover there exists R > such that fl CC Bn, and < faix) < Ce~^l'^ 
whenever |a;| > R. 
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Now for every S > 0, it follows from Theorem 2 and Remark 4 in j4j that 
there exists a sequence C/f £ H^{Br+i]<C) such that y^^jU^ vg weakly in 

L^^^Br+i), and 

(4-14) [lldUff + ^(iC/fP - 1)^) - + \dvs\) 

weakly as measures in -B/j+i. Let e' := cqc, and let := fexU^i, where x S 
C^(5/f+i) is a function such that x = 1 on Bji and jdxl < C- Also, set 
on M'' \ Then as in (j4.7p . and recalling that $e = | loge|<f>; we have 

1 . ^,5. 



The second term on the right-hand side converges to /jja XP^ ' '^s = /ij3 P'f ■ as 
e — )■ 0. The proof of this statement is like that of (|4.12p . but easier. We break the 
other term into two pieces. The first is 

f2 



^X/?(5i«.'.r + j^(it'5i=-i)') 



|loge 

(1 + 0(1)) 



- Iloge' 

Then, since supp(p) C CC Br, it follows from (|4.14p and the uniform convergence 
— > p that 

limsup— ^i?,(xC/f,;Bi?) < / + |dt,5|). 
e-i-o |ioge| J 

And from properties of x and exponential smallness of outside of iJ^j, it easily 
follows that Hf^{xU^,;M.^ \ Br) — > as e — > 0. By combining the above inequalities, 
we find that 

1 

e^o"' |loge|2' 

Note also that 



limsup —^{g,{ul) - g,{f,)) < g{vs) < g{v) + cs. 



(1 + o(l))(/, x)'^^ ^ P«5 weakly in 



|loge| ■ |loge' 

Conclusion (iii) now follows by setting := for (5(e) converging to sufficiently 
slowly. □ 

We conclude this section with the proof of the approximation lemma used above. 

Proof of Lemma \lS\ We introduce some auxiliary functions. First, for r G (0,1], 
let 

{x G : dist(a;, n) < r}, fl^r := {x G O : dist(a:, dfl) > r}. 

Next, for sufficiently small cr > 0, let ^'o- : il^ ^ ^ he the diffeomorphism 
given by 



X if X G 
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where v{x) is the outer unit normal to dil at the point of dfl closest to x, and d{x) 
is the signed distance (positive outside ft, negative in il) to dfl. Note that 

(4.15) \\D^a-I\\oo<CV^, \\D-f-^ -I\\oc<C^ 

where / denotes the identity matrix. In addition, p(5'o-(x)) > p{x) for all x, when- 
ever a is sufficiently small, since Dp[y) = —c{y)v{y) for y S dVL, with c{y) > c > 
for all y, by 

Next, let Xa ^ C^iM?) be a nonnegative function such that x<t = 1 in ^a/A 
and Xcr has compact support in Q,^,^/^^. Finally, for r G (0, 1] let rj^ be a smooth 
nonnegative even mollifier with supported in B{Q,t) with J 77^ = 1. 

We define vs ■= rjr * {xa ■ "^fv), where t, cr will be fixed below. Here ^fv 
denotes the pullback of v by 5*^, which is a one- form on ficr, and the product 
X(j ■ ^'ct w is understood to equal zero on R"^ \ fJ^. 

Note that Xo- • '^fv is integrable on M.^, so that the convolution in the definition 
of vs makes sense. Indeed, (|4.15p implies that |'I'^i;(a;)| < {1 + C y^)\v{'if a{x))\ for 
all X, so that by a change of variables. 



/ \Xa-^*v\< [ Mv\<C [ 

JR3 "'n.1„/4 Ji'c 



\v\ < 00. 



The final estimate follows from the v £ L^, a.s well as the fact that p is bounded 
away from in '^ai^sa/i): since this set is compactly contained in ft. 
We will take r < (t/4, so that vs = rjr * ^fv in f2. Then 

p\vs~v\'^<2 / p|?7^ * ^-fu - ^-fup 4- 2 / p\^*v-v\'^. 



The definition of ^'o- implies that '^fv — v = if dist(j;,R3\f]) > cr, and I**?;-?;! < 
C\v\, so that the second term on the right-hand side tends to as cr —)■ 0, by the 
dominated convergence theorem, and can be made less than 5/2 by choosing cr 
appropriately. Then we can clearly make the first term on the right-hand side less 
than (5/2 by taking r smaller if necessary. 

To estimate / p\dvs\, we consider the action of dvs on some G C^(A^il). 
Exactly as in p.l4p . we can rewrite 



p4>-dvs= vd* *i^^^)*{Tjr*ip ■ 

Since d* ★ {'i'"^)'^ {rjr * (p * 4>)) is smooth, and thus continuous, 

p(t>-dvs = / dv ■ ■k{'i/-'^)* {ijr * {p -k (f))) 

where the right-hand side indicates the integral of the continuous function *(^'^^)'^(- • • ) 
with respect to the measure dv. It follows that 



p 4> ■ dvs < sup 



i(*;l)#(7y,*(p*<^))(x) 



p\dv\ 
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And for x G ^'<t(^1), since | * 4>{x)\ = \(j){x)\ for all x, 

p 



- p{x) ivr*\p*m^^\^)) 



(4.16) <\mJA±^( sup p). 

And writing y := ^^^(cc), if a is small enough then Dp ■ v < —c < in the set 
where Vfo- is not the identity, so it follows from the mean value theorem and the 
definition of that p{'^cr{y)) > piy) + cV^{d{y) + \/o' — (T)~^, where (•••)"'' denotes 
the positive part. Thus 

^ f \ . p{y) + CT 

' sup p < 



p(*a(y)) \B^(^y) J ~ p{y) + c^/^{d{y) + - ct)+ " 

We insist that a < 1/16 (in addition to other smallncss conditions), so that in 
Q,\fl-a-j we have the inequality d{y) + ^/a — cr > sja — 2a > \\fa- In this set, then, 

p[y) + Ct p(y) + Ct ^ rr. ■ . 

< ^\"' < 1 for T sufhciently small. 



p{y) + c^{d{y) + ^-a)+ " p{y) + (c/2)a 

By taking r small enough, wc can make Ct/ p{y) as small as we like in the set fi-o-, 
where p > ca. Thus by taking r still smaller, if necessary, we can guarantee that 
the right-hand side of ()4.16p is bounded by (1 + (5)||0j|oo- Inserting this into the 
above estimates, we conclude that 



p(l> ■ dvs < \\(l)\\oo{l + 5) / p\dv\ 
Jn Jn 

for aU (f) e C^{A'^n), and hence that vs satisfies (|4?T3l) . □ 

4.2. a dual problem and critical forcing. In this section we give the proof of 
Theorem [5l Wc will use notation introduced in Section FOl 



Proof of Theorem\^ 1. Wc first formulate a dual problem. Let 

Io{w):=- I pw<i> + l I p\dw\ = -(w,$)p + i / p\dw\ 
Jn ^ Jn ^ Jn 

so that g{v) = Io{v) + ^\\v\\l. Then 

/o(w) := sup {{v,w)p- Io{w)) = sup ( (w + $, ■u;)p - i [ p\dw\ 



It is clear that if [v + $, 'w)p ~ \ p\dw\ > for any w, then the supremum on the 
right-hand side above is unbounded, so we conclude that 



iiv + ^eN (defined below) 
if not. 



loiv)-- 
for 

(4.17) iV := |e e LliA'n) : {^,w)p < \ j^^P\M for all w G LI{K^^)^ 
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Then it follows from basic facts about duality, see Lemma |6l that the unique min- 
imize! vo of Q in is also the unique minimizer of 

(4.18) Q^v) := Ia{~v) + -\\v\\^ \^ 

i I +00 it not. 

2. We next rewrite the dual problem. It is immediate from the definition ()4.17p 
of N that N C (kero?)^. Hence, writing Pp for L^-orthogonal projection onto 
(ker c?)p, it follows that 

<^ -V e N if and only if Pp* = PpV and P^<^ - P^v € N 

In particular, 

lMl=l\\PrMl + \\\P^vr ii<f-veN 

so that minimizing the norm of v, subject to the constraint $ — w G iV, is 
equivalent to minimizing the norm of P^v, subject to the constraint P^^ — 
P^v e N. 

Now recall from the description (|1.24p of (kcrd)^ = Imagc(P;j'-) that every el- 
ement of (ker(i)p can be written in the form ^-J^ for some d* {3 ^ H]q{A^^l). In 
particular, if we write P^<^ = ^-j^ and PpV[) = then 

vq = Pp^ -\ — where /3o minimizes 



(419) 



^ P {. P 



As usual, we understand ||^^||p to equal +00 if does not belong to Lp(A^f2). 

We now rewrite the constraint by noting that for any smooth w G C°°{A^d) and 
for ^ e Hlf{A^n) such that e L^A^n), 

f^^^~^\ w)p = [ d*{P^-P)-w = / {P^-P)-dw. 



Here the boundary terms vanish due to the fact that (/3$)Ar = /3jv = 0. Then the 
definition (|4.17l) of and facts about density of smooth functions established in 
Lemma [13] imply that 



eN ^ i^-^ — ^,w) < - f p\diu\ for all w e C^{A^n) 
P 2 

[P^ - P) ■ dw <]- I p\dw\ for aU w e C°^{A^Q.) 



P 



(4.20) ^ 11/3* - /3||p* < ^ 



2 

2 



where we recall the definition 

||7||p* := sup{ / "f-dw: w€C°^{A^n), I p\dw\<l}. 



Now by combining (|4.19p and (|4.20p , we obtain the characterization of wo appearing 
in Theorem 3, see ([L^ . p:^ . (HHT]). 
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Observe further that, by stationarity of ()1.22l) with respect to variations 1 1— )■ e^vo 
around t = we obtain 

(4.21) i / p\dvo\ + [ pvn- (uo - $) = . 

^ Jo, Jn 

Recahing that $ = Pp<i> + we have p{vo - $) = (i*/3o - Inserting in 

(|4.2H) yields (|1.28l) after integration by parts. 

3. It remains to check that dv^ = if and only if ||/3$||p* < i. 

The global minimizer of the functional ^ i-^ ||| — |!p in Hj^iA^fl) is attained by 
(3 — 0; and this satisfies the constraint (j4.20p if and only if ||/3$||p* < i. 

Thus if 1 1/3$ lip* < i. then uq = Pp$ G (kerfi)p, and in this case clearly dv^ — 0. 

On the other hand, if ||/3$||p* > i, then vq — Pp^ = is a (nonzero) element 
of (kcrd)^, and hence in this case ^ d{vQ — Pp^) = dv^. 

□ 

5. FURTHER REMARKS 

5.1. symmetry reduction. In the presence of rotational symmetry, the function- 
als we study in this paper reduce to simpler 2-dimensional models. We discuss 
this first for the functional Q, defined in (|1.22p . arising in case of Bose-Einstein 
condensates. 

Lemma 14. Consider the functional Q{v) = /q P w • $ + ||(iw|^ , and as- 
sume that there exist some C [0, oo) x R and some p : — > (0, oo) such that 
ft — {{r cos a, r sin a, z) : (r, z) G a G M}, 
p(r cos a, r sina, z) = p(r, z) Va G K. 
Assume moreover that there exists some : f2 — > M such that 

$(r cos a, r sin a, z) = (/)(r, z)d9 for all a. 

Then the unique minimizer vq of Q is given in cylindrical coordinates by vq — 
WQ{r, z)d9, where Wq minimizes the functional 

(5.1) g'-'^^iw) ■=\j^^P (|Vw| + drdz 

in the space of functions w : — > M such that j^^ juP' dr dz < oo. 

We set g^'^'^{w) = +0O if dw is not a Radon measure in fl or if rp is not \dw\- 
integrable. 

As noted in the introduction, Q^''^'^ is exactly a (weighted) version of a functional 
that has been studied in the context of image dcnoising, see for example [23l [7] . 

Proof. 1. Let Ra : K'^ — > M'' denote rotation by an angle a around the axis. 
Equivalently, in cylindrical coordinates, Ra is the map (r, 6, z) i— > (r, 9 + a,z). Then 
our assumptions imply that 

Ra{fl) poR^ = p. = $ 

for all a. It easily follows that GlR'f^v) = Q{v) for all v and a. By uniqueness of 
the minimizer vq oi Q, which follows from strict convexity, we conclude that 

vo e Ll{K^n) [v e L^iA^n) : R*v ^ v for all a} . 
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It is then immediate that vq minimizes G in Lp(A^f2). 

2. Any v G L^{A^n) can be written in polar coordinates as 

(5.2) V = v^{r,z)de + v'\r,z)dr + v''{r,z)dz. 
We claim that for any such v, 

(5.3) g{v) > Giv'^de) - 2ng'-^''{v') - C($). 

Clearly, this together with Step 1 implies the conclusion of the lemma. To prove 
(|5.3|) . note that if v is smooth and has the form (|5.2p . then 

\dv\ = laydi^ Ade + Oydz Ad9 + {drv' - d,v'')dr A dz\ > \d{v'^d0)\. 

In the general case, the same conclusion follows from the density of smooth functions 
in w G Lp(A^ri) in the sense of Lemma [T51 It is also clear that if v G Lp(A^f7) and 
<i> = (j){r, z)d6, then 

pointwise, so that 

1^ ^(M! _ , . $) = y^J (|. - $p - \^-) > p(^ - iv^de) ■ $). 

Combining these estimates, we conclude that G{v) > Q{v^d9) for all v G L'^{A^il). 
Finally, the identity 

g(w^d0) = 2^g''^'^(/)-C($), C($):= / 

is clear if is smooth, and in the general case can be verified either by a density 
argument similar to the one given above, or by directly relating the definitions of the 
total variation measures associated with d{v^d6) in O and Vu® in (l, respectively. 

□ 

Remark 8. As in Theorem [SJ one can use duality to rewrite the problem of mini- 
mizing Q'^'^'^ as a constrained variational problem. For example, one can verify that 
vq minimizes Q^'^'^ if and only if it minimizes the functional 

(5.4) w 1-^ -w^ drdz 



subject to the constraint 

(5.5) I ^U-wY dr dz <l [ p\VC\ for ah C € C°°(i)), 

analogous to (j4.18p . One could also reformulate this as a problem of minimizing 
a weighted Dirichlet energy of a 1-form on n with a nonlocal constraint like that 
of (|1.26p . but in this setting this seems to us less natural, since the formulation in 
terms of functions rather than 1-forms seems simpler. 

Remark 9. For velocity field represented by the 1-form v = voir, z)dO, the associated 
vorticity 2-form is dv = drW dr Ad9 + dzW dz Ad9. The vorticity vector field, that is, 
the vector field dual to dv, is then ■^{drW iz — dzW ir), where iz and ir denote unit 
vectors in the (upward) vertical and (outward) radial directions respectively. It is 
natural to interpret integral curves of this vector field as "vortex curves" . Since the 
vorticity vector field has no eg component and is always tangent to level surfaces 
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of w, we conclude that, formally, vortex curves have the form "0 = constant, w = 
constant" (at least for regular values of w). 

Thus in the reduced 2d model, we interpret level sets of a minimizer wq, or more 
precisely sets of the form (?{(r, z) : woir, z) > t}, as representing vortex curves. 

For similar reasons, one should think of the "vorticity measure " as being given 
by V^wo, rather than Vwq- 

Similarly, we have 

Lemma 15. Consider the junctional 



and assume that there exist some C [0,oo) x M and : — > M such that 



Then the unique minimizer (vq^Aq) of J- is given in cylindrical coordinates by 
{vo,Ao) ~ {wo(r, z)d6,bo{r, z)dd) , where (u)o,&o) minimizes the functional 



in the space of {w,b) : — > for which J^^'''^{w,b) is well-defined and finite, 
where := {(r, z) : r > 0}. 

We omit the proof, which is extremely similar to that of Lemma [TH 

5.2. contact curves and vortex curves. It is interesting to ask whether one can 
define a useful analog of the "contact set" (as normally defined for classical obsta- 
cle problems) for the variational problems with nonlocal constraints formulated in 
Theorems [5] and We address this question first for Bose- Einstein condensates 
in the presence of rotational symmetry, as discussed immediately above. Thus, we 
assume that wq : Q ^ M minimizes the functional (j5.4p subject to the constraint 
(|5.5|) . An approximation argument starting from (|5.5p shows that if is a set of 
locally finite perimeter in f2, then 



We say that dE is a contact curve if equality holds in the above (where dE should 
be understood as the 1-dimensional set that carries |Vl£;|). 

Lemma 16. For a.e. t, d{wQ > t} is a contact curve. 

As argued in Remark |9l it is natural to interpret d{wQ > i} as a "vortex curve", 
so the Lemma states, heuristically, that every vortex curve for is also a contact 
curve. 

Proof. By using rotational symmetry to reduce (|1.28|) to the (r, z) variables, or by 
using the fact that = we find that 




= {(r cos a, r sin a, z) : (r, z) G 51, a G E}, 
Aexi^ cos a, r sin a, z) = 0(r, z)d9 for all a. 



(5.6) 




(5.7) 





pIVwoI + / -{wq — 0)wo dr dz ~ 0. 
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Using the coarea formula, we rewrite this as 

(5.8) J p|Vl{^„>t}| + J ^{wo-(j))l{^„>t}drdz^dt^O. 

It foUows from (15.71) that 



^ J P|Vl{^o>t}| + J ^{wo ~ (f>)l{nio>t} dr dz > 
for every t, and then (|5.8p imphes that in fact equahty holds for a.e. t. □ 

It is almost certainly not true that every contact curve for the minimizer wq is 
also a vortex curve, in the generality that we consider here, due to the possibility 
of degenerate (nonlocal) obstacles, as in the classical obstacle problem. One might 
hope, however, that the vortex curves and contact curves coincide under reason- 
able physical assumptions (for example, $ = r^d6, corresponding to rotation of a 
condensate around the z axis, probably also with some conditions on p.) 

The situation is more complicated for Bose-Einstein condensates in a general 
domain 51 C M'^ without rotational symmetry, since in this case the analogs of 
vortex curves and contact curves may not in fact be curves and do not in general 
admit a very easy concrete characterization. Abstractly, they may be described as 
follows: if we write Z to denote the closure (in the sense of distributions) of 



{da:ae L^{A^n), [ p\da\ < 1}, 
Jn 



n 

then one can think of the set extr Z of extreme points of (the convex set) Z as anal- 
ogous to the objects — distributional boundaries of sets of finite weighted perimeter 
— used above to describe vortex and contact curves. Indeed, by arguments exactly 
like those of Remark 3 of [37] , general convexity considerations and a bit of func- 
tional analysis imply that extr 2^ is a nonempty Borel subset of a metric space, and 
for any T in the vector space generated by Z (that is, the space U\yoXZ), there is 
a measure fix on extr Z such that 

(5.9) T= Lo d^iriuj) 

J cxtr Z 

and in addition 



(5.10) / pd\T\ = / / pd\u\ dfiTico). 
Jn J cxtr z \Jn / 

We remark that in the closely related situation of divergence-free vector fields on 

M", a concrete characterization of elements of the analog of extrZ as "elementary 

solenoids" is established in [?7] . 

With this notation, an analog of Lemma [16] is 

Lemma 17. Let f3o be the minimizer of the constrained variational problem (|1.25p . 
(|1.26|) ). so that vq = + ^^—^ is the minimizer ofQ{-). Then 

(5.11) / il3^ - I3q) ■ dio < I f pd\uj\. 

Jn ^ Jn 

for every w £ Z. We say that lo G extrZ is a " generalized contact curve" if the 

above holds with equality. 

Furthermore, let pdvo denote a measure on extr Z satisfying (|5.9|) , (|5.10|) ( with 

T replaced by dvo)- Then pdvo o,.e. lo is a generalized contact curve. 



VORTEX DENSITY MODELS 



33 



The proof is exactly like that of Lemma [161 except that (|5.9p , (|5.10p are substi- 
tuted for the coarea formula. Then (|5.1ip follows immediately from the fact that 
(3q satisfies (|1.26l) . and the last assertion us a consequence of p.28|) . 

A version of Lemma [17] could be formulated for the functional J- arising in the 
description of superconductivity and the associated contained variational problem 
described in Theorem [21 using a measurable decomposition ()5.9p , (|5.10|) of the 
vorticity dvo to deduce from (|1.12|) a precise form of the assertion that every (gen- 
eralized) vortex curve is a (generalized) contact curve. 

It would presumably be rather easy to adapt results of [57] to the closely re- 
lated situations considered here, to obtain concrete descriptions of extr Z, or the 
corresponding objects relevant for superconductivity, although we are not sure that 
this would add much insight. It would also be interesting to know whether, if we 
consider the model case of uniform rotation about the z axis (for Bose-Einstein) or 
a constant applied magnetic field (for Ginzburg-Landau), the complexities sketched 
above do not in fact occur, and the vortex curves and contact curves for minimizers 
can in fact be identified with curves of finite length; this seems likely to us to be 
the case. 
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